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Abstract We discuss how different approacties to seiecting member stars in kinematic samptes 
of dwarf spheroidal galaxies affect the estimates of their mass and anisotropy of stellar orbits. 
We demonstrate that the selection of members is an additional source of error compared to 
the usual uncertainties due to the sampling of velocity moments. As an example we use the 
kinematic data set for 202 stars in the Fornax dwarf galaxy for which we model the velocity 
dispersion profile and estimate the mass-to-light ratio and anisotropy assuming that mass 
follows light. We also show that stronger constraints on these parameters can be obtained if 
kurtosis of the velocity distribution is included in the analysis. Using the Besancon model of 
the Milky Way we demonstrate that the majority of contamination in Fornax probably comes 
from the Milky Way stars. 



1 Which stars are members? 

Dwarf spheroidal (dSph) galaxies of the Local Group are believed to be the 
most dark matter dominated galaxies. One of the major problems in determin- 
ing their masses and mass-to-light ratios lies in the selection of members which 
we can trust to trace the dynamics of the galaxies reliably. Due to their prox- 
imity to the Milky Way (MW) the galaxies are strongly affected by tidal forces, 
their stars are stripped and form tidal tails. Most of the stars in the tails are 
not bound to the dwarf and including them in kinematic modelling may over- 
estimate the inferred mass-to-light ratios dramatically. Another possible source 
of contamination are the stellar populations of the MW which can however be 
dealt with using detailed photometric studies (Majewski et al. 2000). On the 
other hand, the tidally stripped stars will have the same photometric properties 
as the population of the dwarf. It is therefore essential to reject them via other 
methods. 

We illustrate a number of methods to deal with such interlopers using a 
recently published kinematic sample of 202 stars in the Fornax dwarf (Walker 
et al. 2006). The sample may contain both tidally stripped stars and those from 
the MW since the selection has been done only by choosing stars along the red 
giant branch and no colour-colour information was used. Figure 1 shows the 
line-of-sight heliocentric velocities of the stars as a function of their projected 
distance from the dwarf centre. In each panel the filled circles correspond to 
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Figure 1: The kinematic sample of 202 stars in tlie Fornax dwarf after applica- 
tion of different schemes for the removal of unbound stars (open circles). 



the stars adopted as members and open circles are interlopers. We present the 

results of four interloper rejection methods described in detail in Wojtak et al. 
(2007): the method based on finding the maximum velocity available to a star 
proposed by den Hartog & Katgert (1996) (HK), the method based on the ratio 
of two mass estimators proposed by Perea et al. (1990), the method based on 
iterative rejection of stars outside ±3(T(i?) (where a{R) is the velocity dispersion 
in different radial bins) and a simple constant cut-off in velocity which retains 
a similar number of stars as the other methods. In the upper left corner of each 
panel we give the number of members and rejected stars. 

2 Modelling of velocity moments 

The velocity dispersion profiles calculated from the different samples are shown 
in Figure 2. As we can see, in the two lower panels the velocity dispersion 
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Figure 2: Velocity dispersion profiles for the different selections of member stars 
in Figure 1. 

profiles show a secondary increase at larger R signifying incomplete removal 
of unbound stars. The velocity dispersion profiles have been fitted with the 
solutions of the Jeans equation (see Lokas 2002, Lokas et al. 2005) assuming 
that mass follows light. The parameters of the light distribution were adopted 
from Irwin & Hatzidimitriou (1995). For each sample we estimated two con- 
stant parameters, the mass-to-light ratio M/Ly and the anisotropy parameter 
P which describes the type of orbits the stars follow. The results in terms of the 
best-fitting parameters (dots) and la, 2a and 3a constraints on the parameters 
(contours) are presented in Figure 3. As we can see, going from the upper left 
to the lower right panel, more and more negative f3 values (corresponding to 
more tangential orbits) are preferred. This is the consequence of the increasing 
velocity dispersion profiles. If the assumption of mass following light were re- 
laxed this type of behaviour could also lead us to overestimate the mass (since 
extended mass distribution can also fit the increasing aios profile). 
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Figure 3: The Icr, 2cr and 3(T constraints on the parameters (3 and MjLy from 
fitting the velocity dispersion profiles of Figure 2 with the assumption that mass 
follows light and /3 =const. Dots mark the best fits. 



The method of HK has been shown by Wojtak et al (2007) to be the most 
effective in removing unbound particles from samples generated from simulated 
dark matter haloes. It has also been shown to perform well on the stellar samples 
drawn from simulated two-component models of dSph galaxies by Klimentowski 
et al. (2007). We can therefore conclude that the results obtained for the HK 
sample are the most reliable. One must however keep in mind when modelling 
the kinematics of dSph galaxies that in addition to the uncertainties due to 
sampling errors of velocity moments there is always additional uncertainty due 
to the method by which the member stars were selected. Since the simple cut-off 
in velocity is still the most common method used for interloper rejection many 
results are biased towards tangential orbits or extended mass distributions. For 
example, Walker et al. (2006) using the same data concluded that if mass follows 
light in Fornax then isotropic orbits (/3 = 0) arc excluded while we show that 
the data are perfectly consistent with /3 = if the unbound stars are properly 
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Figure 4: The kurtosis of the Une-of-sight velocity distribution in Fornax (left 
panel) and the Icr, 2a and 3<t confidence regions on the fitted parameters from 
the joint fitting of dispersion and kurtosis (right panel). 

removed from the sample. 

Once the interlopers are reliably removed one can also consider with con- 
fidence the higher moments of the velocity distribution. Figure 4 (left panel) 
shows the profile of the variable k = (log kios)^^^° where Kios(-R) = vf'^g{R) / af^^g{R) 
is the kurtosis of the line-of-sight velocity distribution. The kurtosis was cal- 
culated for the most reliable sample generated by the HK method of interloper 
rejection. The fourth velocity moment vf^^{R) is governed by the higher-order 
Jeans equation (see Lokas 2002, Lokas et al. 2005). The solid line in the 
Figure is the best-fitting solution when the dispersion and kurtosis are fit- 
ted simultaneously. In the right panel we plot the la, 2a and 'Sa confidence 
regions on the fitted parameters. Comparing with Figure 3 we see that the 
anisotropy parameter is now much more constrained, at 68% confidence we get 
/3 = -0.03t° j3, M/Lv = ll-4li;? Mq/Lq with x^/N = 10.9/10 while when fit- 
ting only the dispersion profile (for the same HK sample) we had /? = — O.IT^q 
M/Lv = ll-3t?:8 M0/L0 with x'^/N = 3.4/4. We sec that when kurtosis is 
added the best fitting anisotropy is even closer to /3 = 0. 

3 The origin of unbound stars 

As discussed in detail by Klimentowski et al. (2007) using an N-body simula- 
tion of a dSph galaxy, the contribution from the tidally stripped stars to the 
kinematic samples depends strongly on the orientation of the tidal tails with 
respect to the line of sight. If the tails are along the line of sight then the con- 
tamination is strongest. This kind of contamination should also be symmetric 
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Figure 5: Expected contamination of the Fornax kinematic sample by MW stars. 

in velocity: we should sec similar amount of unbound stars at velocities above 
and below the mean velocity of the dwarf galaxy. As Figure 1 demonstrates, 
this is not the case: there are much more stars rejected with velocities below 
the systemic velocity of Fornax. There is also evidence for the visibility of tidal 
tails in Fornax from the photometric study of Coleman et al. (2005) which 
means that they cannot be oriented exactly along the line of sight. These two 
facts lead to the conclusion that the majority of the contaminating stars in the 
sample come from the MW stellar population. 

We have verified this conjecture by referring to the Besancon model of the 
MW (Robin et al. 2003). The results are illustrated in Figure 5 where we plot 
the histograms showing the expected number of stars with a given velocity. The 
open histogram refers to all stars expected to be seen in the direction of Fornax, 
the half-filled one to stars along the red giant branch of Fornax and the filled 
one is for stars expected to be seen in a small kinematic sample of 202 stars. 
As expected, the maximum of the distribution falls around the velocity of 20 
km which confirms that most of the contamination in Fornax indeed comes 
from the MW stars. 
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